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ABSTRACT 



Arc tracking is a phenomenon that occur when electric arc occurs 
between two conducting wires whose insulation has broken down. It 
is a major problem in aerospace industry, because of the insulation 
breakdown of Kapton, a popular aerospace insulators. 

It has been shown that ordinary circuit protecting devices, do 
not protect against arc tracking. - This project is primarily 
concerned with a design of a smart circuit to detect arc tracking, 
using the salient features or signature of arc tracking obtained on 
previous studies. 

The smart circuit consists of an inductive pick-up coil, analog 
to digital converter, TMS320C30 processor and relay assembly. The 
processor of the EVM use Welch method to obtain the power spectral 
estimate of the digitalized signal. Then computes the mean value of 
power spectral estimates and the standard deviation of the time 
domain data. These values are to be compared to the previously 
obtained data. When the mean value is less than 30 and the 
standard deviation is greater than 2, the processor is assumed to 
have detected a non-arc tracking signal. Any other conditions 
implies arc tracking. The processor can then activates an 
interrupt. The latter can be used to energize a relay, which will 
open the circuit in cases of arc tracking signal detection. 
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CHAPTER 1 



INTRODUCTION 



m 



Ll ARC TRACKING PHENOMENON 

Aro tracking is a phenomenon that occurs when electric arc 
appears between conducting wires whose insulation has broken down. 
During arc tracking carbon residue is formed across the insulator 
surface due to excessive heat [1]. 

The two basic types of arcing are line fault arcing and contact 
arcing. An arc due to line fault results from either a line-to- 
line short or a line-to-ground short. When fault of this type 
occur, the arc tracking circuit will draw current up to or beyond 
the capacity of the circuit, giving rise to an excessive heat which 
might burn out the insulation and accumulate carbon residue across 
the insulator surface, creating a short circuit. Contact arcing 
results from a high resistance connection in series with a load. 
This may occur due to loose connections, oxidized contacts, or 
foreign non conducting material interfering with conducting path. 
Contact arcing is also commonly caused by use of extension cords of 
insufficient current-carrying capacity. 

Arcing can be gradual with very low characteristics that are 
not detected by ordinary protective devices, which operates within 
a range 15 amps to 3 amps . When an arc conducts current it 
produces high frequency noise on the power line. Tests performed 
by Lockheed corporation indicated that: 
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1.. Circuit protection devices as a rule do not protect 
circuits against arc tracking. However, no arc tracking was 
observed when circuit protection was provided by a 3 amps rated 
remote power control. 

2. Re-closing of the tripped circuit protection devices re- 
initiate the arc tracking phenomenon. 

3, Arc tracking of wire^ in a bmidle cause damage to adjacent 
wiring. 

Figure 1.1 provides a summary on series of tests depicting distance 
of arc track for each protection devices. |f| 

1.2 PRELIMINARY STUDIES 



(FFT) was used. The parameters obtained from arc tracking studies 
are shown in Table 1.1 and 1.2 

The author found that [3,4] , of all of all parameters shown in 
Table 1.1 and 1.2, only the mean value of the power spectral 
density, the standard deviation of the power spectral density of 
the time domain data, and the percentage arc tracking occurred can 
be used to differentiate arc tracking from non-arc tracking 

2 
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In a study conducted by August Kowero, John O. Attia and 
Dwight Janoff [3], Kapton insulated signal cables which were made 
bare at the end, were brought together in a container containing 
graphite and copper powder. Power supply of 20 V, 3 V, and 42 V i|| 
were applied and current was varied. The Welch method for direct 



iiii 



i||il8 

Mi 

computation of power spectral density using Fast fourier transform ||| 



m 

ifl 



tfi'^ 



1. 11 






00^58 



Case 1 :06-cv-00584-RMC Document 6-4 Filed 06/07/2006 Page 9 of 25 



as 

CO 



20 A CB ' 
20 A FUSE 

20 A RPC 
15 A FUSE 
10ACB \\ 
10AFUSE — J 



I ♦ 

4 3* (#12) 



I 7*(#12.#10.#8) 



HH 3* (#12) 



+ 



-I 3' (#16) 



10 A RPC 

7 AGO 

8 A FUSE 

7.5 A RPC 

5ACB 

5 A FUSE 
5 A RPC 

3ACB 
3 A FUSE 

3 A RPC 



«| 3' (#20) 



-I 3* (#22) 



-I 7* (#22.#24) 



-16* (#22.#24) 
1 3* (#26) 



6* («24.#26) 



6* (#24.#26) 
h— 



+ 



+ 



+ 



r 



2" 3" 4" 

INCHES OF ARC TRACK 



-I 15" (#16.#20) 
-I 6* (#20) 



— M 3* (#22) 
\ 3* (#22) 



-110* (#22,#24) 



I I 



+ 



5" 



6" OR MORE 



I .MINMAXaURNLENGTHE 

• .AVERAGE OFALLTESTS 

• -TOTAL NUMBER OF TESTS 
l»). WIRE SEES TESTED 



FIGURE 1.1 SUMMARY OF 28VDC ARC TRACK TEST WITH PROI^ER 
SIZED CIRCUIT PROTECTION S 



Case 1 :06-cv-00584-RMC Document 6-4 Filed 06/07/2006 Page 1 of 25 



TABLE 1.1 



Parameters Obtamed from Arc Tracking Data 



standard deviation of time domain arc tracking data 
peak value of power spectral density 
mean value of power spectral density 
standard deviation of power spectral density 



SD - 

MP - 

MVP 

SDP - 

MFRQ - Mean frequency content 
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signals. The percentage .of time are tracking occurred was found to 
be not useful for real time detection of arc tracking. Using the 
two parameters in fable 1.1, standard deviation of the time domain 
data and the mean value of power spectral density Figure 1.2 was 
drawn. 

From Table 1.1 and Figure 1.2, it can be concluded that for arc 
tracking to occur, the mean value df^rhe power spectral density 
have to be greater than 3 , while the standard deviation of the 
time domain data can have any value. In addition, a signal can be 
arc tracking provided the mean value of the power spectral density 
is less than 3 and the standard deviation of the time domain 
signal is less than 2.0. 
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1.3 OBJECTIVES 

Because arc tracking is not easily detected by ordinary circuit 
protecting devices, there is a need for a smart circuit that can be 
used to detect arc tracking. The objective of this work is to 
design a smart circuit that can be used to detect arc tracking. 

Chapter 2 introduces the software and tools associated with 
TMS3 2 0C3 EVM. A technique of computing power spectral density are 
discussed. Chapter 3 has the block diagram and the methodology 
used in designing the smart circuit. It summarizes the DSP 
algorithm, and a description of software program written for arc 
tracking detection. The conclusion of this work can be found in 
Chapter 4. In addition, the latter chapter are some comments on 
future work. 
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TABLE 1.2 



Other Parameters of Arc Trackmg Data 

TARC - Percentztge t)f Time Arc Tracking occurred 
BUPT - Percentage of Time for the Current to Build-Up 
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FIGURE 1.2 PLOT OF THE STANDARD DEVIATION VERSUS MEAN 
VALUE OF THE POWER SPECTRAL DENSITY 
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CHAPTER2 



DIGITAL SIGNAL PROCESSOR AND SPECTRAL ESTIMATION 






2.1 ARCHITECTURE OF TI TMS320C30 EVM. 

The DSP application require an analog input and output signal. 
Along the input path is an analog to digital converter, which will 
convert the input signal to digital signal. The digitalized 
signals pass through the signal processor where the analysis of 
whether they contain arc tracking phenomenon or not is done. The 
output of the processor is converted to analog signal by digital to 
analog converter, The output signal is used to trigger a relay for 
the detection of arc tracking parameters. 

The TMS3 2 0C30, shown in figure 2.1, and function block diagram 
on figure 2.2, is a 32-bit digital signal processor, capable of 
performing floating-point, integer and logical operation. It has 
two serial ports, two timers and on-chip direct memory access (DMA) 
controller for I/O. These features are controlled through memory- 
mapped registers, to communicate or to interface with the outside 
world. The serial ports are interrupt-driven, with interrupt 
vectors at specific memory location. 

Along the input path is a filter to suppress frequencies 
greater than one-half the sampling frequency, the analog-to-digital 
converter (ADC) , to convert the input signal to digital signal. 
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FIGURE 2.1 TMS320C30 EVM BLOCK DIAGRAM 
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The maximum level of the input signal for conversion is based on 
the specific ADC, (EVM is +/- 1.5 V ,TI AIB is +/- 10 V ) . 

After processing, the digital-to-analog converter (DAC) 
converts the signal to analog signal. The levels from the DAC are 
produced based on the digital word on its input. Step level signal 
in passed through the low-pass output filter to smoother out the 
steps and reconstruct the equivalent -analog signal based on the 
step levels. 

The methodology utilized in distinguishing arc tracking from 
non-arc tracking, signals is dependent on the presence of salient 
distinguishable arc tracking features [3]. 

The EVM is installed into the PC as an I/O device. The EVM 
data and command register are at different addresses in the PC's 
I/O, with the command register looking like 8 bit I/O, and the data 
transfer register looking like 16 bit I/O. The PC can send, 
receive or command information simply by reading or writing to 
different I/O addresses on its bus even though it is physically the 
same register. 

The EVM analog section consists of TLC 32044 Analog interface 
circuit coupled to an LM 386 low noise audio power amplifier, for 
output drive, and TL 07 2 input amplifier. The ability to provide 
emulation support within target system (embedded emulation) is 
enabled by the SN 74ACT8990 Test bus control and the TMS320C30 
Emulation port. TMS3 2 0C3 can be loaded executable programs via the 
loader without having to get into the debugger {EVM3 0) . The 
command format for the loader is: evmload <-options> filename. out 
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2.2 POWER SPECTRAL ESTIMATION 



One method of power spectral estimation is the Welch method 
[4]. The latter method estimate the spectrum by method of averaging 
periodogram using FFT. Sampling in the time domain is associated 
with the periodicity in frequency domain and vice versa. In Welch 
method the sequence of x(n)* data haS a power spectrum denoted by 
P^(W). Where the Fourier transform of x(n) is discrete, l:he 
fundamental frequency of the signal in the time domain and thus the 
periodicity is maintained as the sum of the harmonically related 
complex exponential sequences, 

The periodogram spectral estimates is obtained by dividing x(n) 
from < n < N - 1 data points into segment K overlapping segments 
of length L. In implementing segments overlapping by L/2 obtaining 

K = {(N - (L/2)} 
The ith segment of the data is defined as 

Xi(n) = x(il/2 + n) W,(n) 



iljl 



where 0<n<i<K-l; 0<i<K-l 
The Wd(n) is an L point data window. 

The M point of the DFT of the window point is given by modifying 
periodogram equation 
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Then it is averaged to produced normalized spectral estiiriate, at 
radian frequency • 

- ^-1 i-i 

S^{2nk/M) ^^J2S^{K) while U ^Y^w'^in) 



for <= k <= -1 

U is necessary for the spectral estimate to be unbiased. Appendix 
A shows a program for power spectral estimation. 

In appendix B, examples of non-arc tracking and arc tracking 
signals in the time domain data and their power spectral density 
are shown. 
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CHAPTER 3 



SMART CIRCUIT DESIGN 



3.1 SMART CIRCUIT BLOCK DIAGRAM 

The use of TMS3 2 0C3 EVM and the implementation of FFT wii.1 
gives the frequency content and time content of the input signal 
x(t) . The design of the smart circuit entails calculations of the 
mean of the frequency domain data, and the standard deviation of 
the time domain data, which are in-turn compared to the previously 
determined values. Based on this logic an interrupt is set 
indicating whether an arc tracking has occurred or not. 

The smart circuit for arc tracking detection is shown in 
figure 3.1. It consists of Texas instrument digital signal 
processor TMS3 2 0C3 Evaluation module (EVM) , an inductive pick-up 
monitor and a relay. 

The pick-up signal from the inductive coil is converted to 
digital signal using an analog-to-digital converter of the 
Evaluation module. The digitalized signal is then analyzed by the 
TMS320C30 processor. If arc tracking features are detected an 
interrupt signal is activated. The interrupt signal can then be 
used to energize a relay that will disable the power supplied to 
sensitive equipment. The intelligence of the circuit lies in the 
processor and the software program of the evaluation module. 
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3.2 DSP ALGORITHM 



The Smart circuit for arc tracking is primarily a decision 
making device working on algorithms whose function is to process 
the input signal to the system and decide what the signal 
properties are pertaining to known specifications and conditions of 
arc tracking signals. * - \ 

The realization of the real-time implementation of the follow 
several basic steps. A 2 kHz sinusoidal signal was from the signal 
generator connected to the EVM. The power spectral data is 
obtained from the fast Fourier transform of the input signal. The 
variance and the standard deviation of the time domain data can be 
calculated as: 
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variance =-j^ g U^ - 



mean) 
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Standard deviation = ^variance 



The mean value is computed as: 
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The program shown in Figure 3.2 is for the real time 
implementation of a 512-points real value input FFT. This program 
call the program shown in Figure 3.3 and Figure 3.4 for the 
computation of the mean value of the data and the standard 
deviation of the time domain data. 

A c program SINEGEN^C was used to generate the twiddle constants 
for a 512-point real valued. input FFT [6j. A total of N/2 or 256, 
points are required (for a complex FFT (5/4) N points are requireST 
for W) . The sinegen.c was run. The 256 twiddle constants for the 
512 points real valued FFT are stored in the file twid512.asm. The 
vecsir .obj file was linked for the reset vector. FFT512r.out was 
run using the EVM by 'EVMLOAD512R.OUT' instruction. 
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/*FFT512R.C-REAL-TIME FFT WITH 512 POINTS. CALLS FFT_RL*/ 

#include <math.h> /*std library fiinc */ ' 

#include "aiccom.c" /*AIC com routines */ 

#defmeN512 /*sizeofFFT */ 

#define M 9 /*number of stages */ 

volatile int index = 0; /*input_output index*/ 

float *IO_buffer, *data, *temp; /*pointers to array buffers*/ 

int AICSEC[4] = {0xl428,0xl,0x4A96,0x67}; /*AIC config data */ 

extern void fft_rl(int, int, float *); /*fft function protype*/ 

volatile int *INTVEC = (volatile int *) 0x00; /*addr of inter vecs */ 

void cJntOSO /*interrupt hafldler function */ 

{ * 

PBASE[0x48] = ((int)(IO_buffer[index])) < < 2; /*output data */ 

IO_buffer[index] = (float)(PBASE[0x4C] < < 16 > > 18); /*input data " 
if (+ +index > = N) index = 0; /*increment index, reset = N */ 

} 

mainO' 

{ 

int loop; 

float real, img; 

INTVEC[5] = (volatile int) c_int05; 

AICSETJO; 

IO_buffer = (float *) calloc(N, sizeof(float)); /*input_output buffer */ 

data = (float *) calloc(N, sizeof(float)); /* fft data buffer */ 

while (1) /* create endless loop */ 

sd (float IO_BUFFER0,resuIt2,int n.float *mv,float*sd,float *var,*IO_OUTPUT)/*cail std dev prog*/ 

fft_rl(N, M, (float *)data); /*call FFT function */ 

data[0] = sqrt(data[0]*data[0])/N; /*magnitude of X(0) */ 
for (loop = 1; loop < N/2; loop++){ /*calculate X(l)..X(N/2-l)*/ 
real = data[loop]; /*real part */ 

img = data[N-Ioop]; /*imaginary part */ 

data[loop] = sqrt(real*reaH-img*img)/N; /*find magnitude */ 

} 

data[N/2] - sqrt(data[N/2]*data[N/2])/N; /*magnitude of X(N/2) */ 
for (loop == N/2 + 1; loop < N; loop++) /*X(N/2 + l).. X(N-1) 



/* declare variable */ 

/* declare variables */ 

/^install inter function "^1 
/*config AIC for interrupt*/ 



V 



data[Ioop] = data[N-loop]; 
while (index); 
temp = data; 
data = IO_buffer; 
IO_b offer = temp; 
IO^buffer[0] - -2048; 



/*use symetry 



/*wait till IO__buffer is empty*/ 
/*temp pointer = > data buffer */ 
/*IO_buffer becomes new data buffer */ 
/*data buffer becomes new lO^buffer */ 
/*frame sync pulse (negative spike) */ 

mv (float IO_BUFFER[J,result,int n,fIoat *mv,*sd,*IO_o"uTPuf); /* call mv.c */ 
prog*/ 

} FIGURE 3.2 PROGRAM TO COMPUTE MEAN AND STANDARD 

} DEVIATION FROM POWER SPECTRAL DENSITY. 
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/* mv.c */ 

/* This routine tests the data and return its' mean value. */ 

/* The input data[l...n] is drawn from IO_BUFFER */ 

/* n is the number of sample data */ 

/* The output data[sd] and data[mv] are store in *mv and then be /* 

/* transfered to the IO_OUTPUT register */ 

/* respectively */ 

/* */ 



^include <math.h> 

void mv (float IO_BUFFERn,int n, result, volatile int *IO_OUTPUT, float *mv,float *sd) 

{ 

volatile int *IO_OUTPUT = (volatile int *) ox804001 

void n_error (char error_testD); 

intj; 

float ep = 0.0, s,p; 



if (n < = 1) /* data test 

n_error("n must be at least 2"); 
s = 0.0; 



"/ 



result = 0; 

forG = l;j <= n; j + +) 

s = s + IO_BUFFER0]; 

*mv = s/n; 

result = *mv; 

*IO-OUTPUT =(int) result; 



} 



FIGURE 3.3 PROGRAM FOR OBTAINING MEAN VALUE 
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/* sd.c ' */ 

/*This routine tests the data and compute standard deviation. The input*/ 

/*data[l...n] is drawn from IO_BUFFER, n is the number of sample data */ 

#include <math.h> 

void sd (float IO_BUFFERn,int n,int result2, volatile int *IO_OUTPUT,float *mv,float 

*sd, float, *var) 

{ 

volatile int *IO_OUTPUT = (volatile int *) o x 804000; 

void n_errcr (char error_testQ); 
'intj; Z' 

float ep = 0.0, s,p; 



if (n < = 1) /* data test */ 

n_error("n must be at least 2"); 
,s = 0.0; 

result2 = 0; 

forO = l;j <=n;j + +){ 

s = s + IO_BUFFER[j]; 
*mv = s/n; 

*adev = (*var)=(*skew)=(*curt)=0.0; /* second pass for ABS */ 
forO =l;j <=n;j + +) 

{ 

*adev + = fabs(s = IO_BUFFER[j] - (*mv)); 

*var + = (p = s*s); 
*skew + = (p *= s ); 
*curt + = (p *= s ); 

} 

*adev / = n; -^ 

*var = (*var - ep*ep/n)/(n-l); 
*sd = sqrt(*var); 
result2= *sd; 
*I0_0UTPUT=result2; 
if (*var) 

{ 

*skew /= (n*(*var)*(*sd)); 

*curt = (*curt)/(n*(*var)*(*var))-3.0; 

} 

else n_error("No skew / kurtosis when variance=0); 

} 
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FIGURE 3.4 PROGRAM TO TEST THE DATA AND CALCULATE THE 
STANDARD DEVIATION OF THE TIME DOMAIN DATA 
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